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We present a theoretical study of dielectrophoretic (DEP) crossover spectrum of two polarizable particles 
under the action of a nonuniform AC electric field. For two approaching particles, the mutual polarization 
interaction yields a change in their respective dipole moments, and hence, in the DEP crossover spectrum. The 
induced polarization effects are captured by the multiple image method. Using spectral representation theory, 
an analytic expression for the DEP force is derived. We find that the mutual polarization effects can change the 
crossover frequency at which the DEP force changes sign. The results are found to be in agreement with recent 
experimental observation and as they go beyond the standard theory, they help to clarify the important question 
of the underlying polarization mechanisms. 
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I. INTRODUCTION 

When a polarizable particle is subjected to an applied elec- 
tric field, a dipole moment is induced into it. The movement 
of colloidal particles in an applied AC electric field is called 
dielectrophoresis [Jl|]. It is typically used for micromanipula- 
tion and separation of biological cellular size particles, and 
it has recently been successfully applied to submicron size 
particles as well. Specific applications include diverse prob- 
lems in medicine, colloidal science and nanotechnology, e.g. 
separation of nanowires 0, viruses [^[j, latex spheres [|], Q], 
DNA [H and leukemic cells [0], as well as lab-on-a-chip de- 
signs [pp. 

The dielectrophoretic (DEP) force exerted on a particle can 
be either attractive or repulsive depending on the polariz- 
ability of the particle in comparison to the medium. For a 
nonuniform AC electric field, the magnitude and the direc- 
tion of the DEP force depends on the frequency, changes in 
surface charge-density and free charges in the vicinity of the 
particle. The frequency at which the DEP force changes its 
sign is called the crossover frequency (/cf)- Analysis of the 
crossover frequency as a function of the host medium conduc- 
tivity can be used to characterize the dielectric properties of 
particles, and is at present the principal method of DEP anal- 
ysis for submicrometer particles [^|, |j. 

In the dilute limit, i.e, when a small volume fraction of 
charged particles are suspended in an aqueous electrolyte so- 
lution, one can focus on the DEP spectrum of an individual 
particle ignoring the mutual interactions between the particles. 
Although the current theory ^ captures some of the essential 
physics in the dilute case, it is not adequate [|[ [l(J [ill [j~2[ fj~3| ] - 
This is due to the fact that even for a single colloidal particle 
in an electrolyte, it is not established which mechanisms con- 



trol its dielectric properties. If the suspension is not dilute, 
the situation becomes even more complicated due to the mu- 
tual interactions between the particles. One should also note 
that particles may aggregate due to the presence of an exter- 
nal field, even when the suspension is at the dilute limit under 
zero field conditions. In this case, the mutual interactions have 
to be included in the description. 

In this article, we present a theoretical study of the DEP 
spectrum of two spherical particles in the presence of a 
nonuniform AC electric field. We use the multiple image 
method [|l4]|, which is able to capture the mutual polarization 
effects. Using the spectral representation theory Jl5[], we de- 
rive an analytic expression for the DEP force and determine 
the crossover frequency. Our theoretical analysis shows that 
the induced mutual polarization interactions plays an impor- 
tant role in DEP spectrum. In a more general framework, our 
results demonstrate the importance of correlation effects. This 
is analogous to the findings in charged systems where phe- 
nomena such as overcharging, or charge inversion (see e.g. 
Refs. flip] , |l7|]), provide spectacular demonstrations of corre- 
lation effects. 

As our starting point, we consider a pair of interacting 
charged colloidal particles dispersed in an electrolyte solution. 
When the two particles approach each other, the mutual po- 
larization interaction between them leads to changes in their 
respective dipole moments [18|, and hence also in the DEP 
spectrum and crossover frequency. We analyze two cases: 1) 
longitudinal field (L), in which the field is parallel to the line 
joining the centers of particles, and 2) transverse field (T) in 
which the field is perpendicular. The former corresponds to 
positive dielectrophoresis where a particle is attracted to re- 
gions of high field and the latter to the opposite case, referred 
to as negative dielectrophoresis. 
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This paper is organized as follows. In Sec. |l| we present 
the formalism and derive analytic expressions for the effec- 
tive dipole factors in spectral representation. In Sec. |Ul| we 
use the analytical results to numerically solve the crossover 
frequency, dispersion strength and DEP spectra under differ- 
ent conditions. This is followed by a discussion of the results 
in Sec. 0. 



II. FORMALISM AND ANALYSIS 

First, we consider a single charged spherical particle sus- 
pended in an electrolyte and subjected to a nonuniform AC 
electric field. The DEP force Fdep acting on the particle is 
then given by | 



Fdep 



i7re 2 I? 3 Re[fe]V|E| 2 , 



(1) 



where D is particle diameter, e 2 the real dielectric constant of 
host medium, E the local RMS electric field, and Re [b] the 
real part of the dipole factor (also called Clausius-Mossotti 
factor) 



ei - e 2 
h + 2e 2 ' 



(2) 



Here, ?i and e 2 are the complex dielectric constants of the par- 
ticle and the host medium, respectively. In order for the two 
above equations to be valid in an AC field, the dielectric con- 
stant must include dependence on the frequency. The complex 
frequency dependent dielectric constant is defined as 



£27r/' 

where e is the real dielectric constant, a denotes conductivity, 
/ the frequency of the external field, and i = 

The conductivity of a particle consists of three components: 
Its bulk conductivity (oibuikX surface effects due to the move- 
ment of charge in the diffuse double layer (conductance kd), 
and the Stern layer conductance (fc s ), i.e., 



Cl = Clbulk 



D 



D 



(3) 



The diffuse double layer conductance kd can be given as [ |20| 

4F; 2 c2 2 S(l + 3A/2 2 ) 



kd = 



R T k 



cosh f — 1 — 1 



where S is the ion diffusion coefficient, z the valency of coun- 
terions, F a the Faraday constant, Rq the molar gas constant, C 
the electrostatic potential at the boundary of the slip plane and 
To the temperature. The reciprocal Debye length k providing 
a measure for screening on the system is given by 



where c is the electrolyte concentration. Parameter A in 
Eq. (Q) describes the electro-osmotic contribution to kd, and it 
is given by 



A = 



RqTq 

F a 



2e 2 
3r?5' 



(6) 



where r\ is the viscosity of medium. In addition, the Stern 
layer conductance k s has the form [ [TT| ] 



2iV 



(7) 



where u is the surface charge density, £ molar conductivity 
for a given electrolyte, and /i r gives the ratio between the ion 
mobility in the Stern layer to that in the medium. 

For a pair of particles at a separation R suspended in an 
electrolyte, we have to consider the multiple image effect. We 
consider two spheres in a medium, and apply a uniform elec- 
tric field Eo — Eqz to the suspension. This induces a dipole 
moment into each of the particles. The dipole moments of par- 
ticles 1 and 2 are given by p 10 andp 20 (=pio = e 2 T -D 3 6/8), 
respectively. 

Next, we include the image effects. The dipole moment pio 
induces an image dipole pu into sphere 2, while p\\ induces 
another image dipole in sphere 1 . As a result, multiple images 
are formed. Similarly, p 2 o induces an image pi\ into colloid 
1 . The formation of multiple images leads to an infinite series 
of image dipoles. 

In the following, we obtain the sum of dipole moments 
inside each particle, and derive the desired expressions for 
dipole factors. We consider two basic cases: 1) longitudinal 
field (L), where the field is parallel to the line joining the cen- 
ters of particles, and 2) transverse field (T), where the field is 
perpendicular to the line joining the centers of particles. Us- 
ing the above ^notation, the effective dipole factors for a pair 
are given by 



b L * = 

oo 

b T * = bj2(-w 



sinh i 



sinh(n + l)a 

sinh a 
sinh(n + l)a 



(8) 



where a is defined via cosh a = R/D. The summations in 
Eqs. (§) include the multiple image effects, the n = term 
giving the dipole factor of an isolated particle. 

We have to derive the analytic expressions for Re[&2] an d 
(4) Re[6J] to resolve the DEP force in Eq. ([!]). To do that, we 
resort to spectral representation theory. It offers the advan- 
tage of being able to separate the material parameters (such as 
dielectric constant and conductivity) from structural informa- 
tion JI^ ] in a natural way. 

Let us begin by defining a complex material parameter s — 
1/(1 — ei/e 2 ). Using this, the dipole factor for a pair takes the 
form 



e 2 i?o7o 



(5) 



F„ 



(9) 



n=l 
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where n is a positive integer, and F n and s n are the n— th 
microstructure parameters of the composite material As 
an example, the dipole factor of an isolated particle in spectral 
representation expression becomes b = F\/(s — s\), where 
Ft = -1/3 and si = 1/3. 

In order to obtain expressions for the dipole factors b* L and 
Kp in Eqs. (m, we introduce the following identity 



1 



sinh x 



Its application into Eqs 
formations: 



bl = 



4m(m + 1) exp[— (1 + 2m)x\. 

yields the following exact trans- 
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(10) 



where the m-th components of the microstructure parameter 
of the composite material are given as 
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2exp[-(l 
{1 + cxp[-(l -+ 



+- 1) sinh 3 aexp[-(2m + I) a], 
+ 2m)a]}, 
2m)a]}. 



To make this approach more tractable, we introduce di- 
mensionless dielectric constant and conductivity [pT|], s = 
1/(1 — £1/62) and £ = 1/(1 — ax/crz), respectively. Now, 
we are able separate the real and imaginary parts of the ar- 
guments in the expressions for b* L and b^ in Eq. (10). The 
argument can be rewritten as 

. Ae m ///„ 



s - s r 
where 



and 
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The analytic expressions for Re[6£] and Re[6^] (Eq. ( 10)) be- 
come 

? {L) 
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Ac 
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1 + P/f, 



W ). (14) 
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Using these, we can obtain the DEP force Fdep which in- 
cludes corrections due to the image effects. The DEP spec- 
trum consists of a series of sub-dispersions with strength Ae m 
and characteristic frequency / mc . In particular, the frequency 
which yields F = 0, namely Re [b*] = 0, is the desired 
crossover frequency /cf- 
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FIG. 1: DEP crossover frequency vs. medium conductivity for an 
isolated particle (solid line) and two particles at different separations. 
L denotes longitudinal field case and T transverse field case. Param- 
eters: C = 0.12 V, 77 = 1.0 x Kr 3 Kg/(ms), u = 0.033 C/m 2 , 
£ = 0.014 Sm 2 /mol, 6l = 2.25e . 



III. NUMERICAL RESULTS 

The above formalism enables us to study the effects due 
to multiple images under different physical conditions and to 
compare the theory to experimental results. In the following, 
we compare the crossover frequency of an isolated particle to 
that of two particles at different separations. We study the ef- 
fects due to multiple images by varying medium conductivity, 
the (^-potential, medium viscosity, surface charge density, real 
dielectric constant of the particle and molar conductivity. Fi- 
nally, we have computed the DEP spectrum and the dispersion 
strength. 

The common parameters used in all numerical computa- 
tions are the following: Temperature To = 293 K, dielectric 
constant of host medium e 2 = 78eo, bulk conductivity of the 
colloidal particle cruuik — 2.8 x 10 -4 S/m, ion diffusion coef- 
ficient 3 = 2.5 x 10~ 9 m 2 /s, the ratio between the ion mobility 
in the Stern layer to that in the medium fi r = 0.35, particle 
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FIG. 2: DEP crossover frequency vs. medium conductivity when the 
particle size is varied. Parameters as in Fig. 



FIG. 3: DEP crossover frequency vs. medium conductivity for dif- 
ferent ^-potentials. Parameters as in Fig. [l]. 



diameter D = 2.16 x 10 7 m, counterion valency z = 1. The 
dielectric constant of vacuum is denoted by eo- 

Figure [j] shows the DEP crossover frequency as a func- 
tion of medium conductivity for an isolated particle and for 
two particles at different separations. In agreement with re- 
cent experiments [gg], we find that a peak in the crossover 
frequency appears at a certain medium conductivity. The ap- 
pearance of a peak is preceded by an increase of /cf upon 
increasing medium conductivity and followed by an abrupt 



drop 1 1 1 



Compared to an isolated particle, the multiple 
image effect leads to a red-shift (blue-shift) in /cf in the lon- 
gitudinal (transverse) field case. Furthermore, for longitudi- 
nal (transverse) field, the stronger the polarization interaction, 
the lower (higher) the crossover frequency. In addition, it is 
worth noting that the effect of the multiple images is the op- 
posite in the longitudinal and transverse cases. As the ratio 
R/D grows, the predicted crossover spectrum approaches to 
that of an isolated particle, i.e., at large separations the multi- 
ple image interaction becomes negligible. 

Motivated by a recent experiment JT^], we analyzed the ef- 
fect of particle size on the crossover frequency by keeping the 



ratio R/D fixed and varying the particle diameter. In agree- 
ment with the experiments, we find that the location of the 
peak is shifted to higher frequencies and higher conductivities 
when the diameter of the particle is reduced, see Fig. ||. 

Figure || displays the effect of the ^-potential. It has been 
experimentally observed by Hughes and Green [|TT| that de- 
creasing the ^-potential may red-shift the DEP crossover fre- 
quency. The system used by them contained many latex beads 
suspended in a solution, and hence the multipolar interaction 
is expected to play a role. Our results are in qualitative agree- 
ment with the above experimental findings. Furthermore, an 
increase in the ^-potential leads to higher /cf in both the lon- 
gitudinal and transverse field cases. Similarly, increasing the 
real part of the dielectric constant leads to an increase in /cf, 
as displayed in Fig. [|. Increasing the viscosity of the medium 
(figure not shown here), however, has exactly the opposite ef- 
fect for both the longitudinal and transverse field cases. 

Figure [| shows the effect of molar conductivity S on 
crossover frequency. For small medium conductivities (here, 
(72 < 10~ 2 S/m), increasing £ leads to an increase in the 
crossover frequency. However, there is a crossover after which 
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FIG. 4: DEP crossover frequency vs. medium conductivity. The real 
part of the dielectric constant is varied. Parameters as in Fig. |l|. 



FIG. 5: The effect of molar conductivity on the DEP crossover fre- 
quency. Parameters as in Fig. |l| 



lower values of £ yield higher /cf- Similar behavior for the 
low surface conductivity regime has been observed in exper- 
iments JTT|], but the authors are not aware of any systematic 
study of the molar conductivity on Jqf- As Fig. g shows, the 
effect is similar for both longitudinal and transverse fields. 

Figure ^ shows the effect of varying the surface charge 
density on the crossover frequency. Variations in the surface 
charge density lead to more pronounced effects in the low fre- 
quency region, but close to the peak the variations differences 
are very small. In addition, the location of the peak is only 
weakly dependent on surface charge density. These results 
are in agreement with the experimental observations of Green 
and Morgan [Q. 

In Fig. 0, we investigate the real part of the dipole factor, 
and thus the DEP force. The figure shows that the effect due 
to multiple image plays an important role at low frequency 
region when the particles separation is not large, whereas its 
effect is smaller in the high frequency region. In the low fre- 
quency region, the DEP force is be enhanced (reduced) due to 
the presence of multiple images for longitudinal (transverse) 
field case. As the particle separation grows, the multiple im- 



age effect becomes negligible as expected. We also studied 
the effect of particle size on the real part of the dipole factor 
and the effect of multiple images increases as the particle size 
decreases, and the effect is stronger in the longitudinal field 
case. 

Finally, in Fig. |l we plot the dispersion strengths (Ae™ ' 

(T) 

and Ae^ ') as a function of the characteristic frequencies 

(L) (T) 

(fmc and fmc), for m = 1 to 100 with different medium con- 
ductivities 02. Here, m is a positive integer, and F m and s m 
are the microstructure parameters of the composite material, 
see Eqs. [Tl|-[l3[ Hence, Ae m and f mc are the m-th dispersion 
strength and characteristic frequency due to the presence of 
multiple images as discussed in Sec. R. 

The advantage of using the spectral representation theory 
is shown in Fig. ||. Based on Fig. 0, it may appear that 
only one dispersion exist. Figure shows, however, that sub- 
dispersions with strength Ae m and characteristic frequency 
f mc co-exist, and most of them lie close to the main disper- 
sion. Thus, the spectral representation theory helps us to gain 
more detailed information about the system and it provides a 
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FIG. 6: The effect of surface charge density on the DEP crossover 
frequency. Parameters as in Fig. [l]. 



FIG. 7: DEP spectrum (the real part of the dipole factor). Parameters 
as in Fig. [l| 



detailed comparison between the longitudinal and transverse 
field cases. 

At a given 02, for the longitudinal (transverse) field case, 
increasing m leads to corresponding sub-dispersions in the 
characteristic frequency due to the presence of multiple im- 
ages. The crossover frequencies f c are 3.49 x 10 6 Hz, 1.0 x 
10 7 Hz and 1.4 x 10 7 Hz (Fig. §). From Fig. § we find that at a 
lower medium conductivity (say, 02 = 1.0 x 10~ 4 S/m), mul- 
tiple images have a stronger effect on the DEP spectrum for 
the longitudinal field case than for the transverse field. This 
is also apparent in Fig. ^ as well. Moreover, for longitudinal 
field case the multiple images play a role in the low frequency 
range (i.e., smaller than f c ). For the transverse field the situa- 
tion is the opposite. At a larger 02 (say, 02 = 5.0 x 10 -3 S/m 
or CT2 = 1.0 x 10 _2 S/m ), the sub-dispersion strengths for 
the two cases have only a minor difference. These observa- 
tions may partly explain the results of Green and Morgan [ pj| ] 
whose data suggests that there exists a dispersion below the 
frequencies predicted by the current theory. The importance 
of these observation lies in the fact that they help to clarify 
the interesting question of which polarization mechanisms are 



present. 



IV. DISCUSSION AND CONCLUSION 



In this study, we have investigated the crossover spectrum 
of two approaching polarizable particles in the presence of 
a nonuniform AC electric field. When the two particles ap- 
proach, the mutual polarization interaction between the par- 
ticles leads to changes in the dipole moments of each of the 
individual particles, and hence in the DEP crossover spectrum. 
This can be interpreted as a correlation effect analogous to the 
ones seen in charged systems 1 16 1. 



For charged particles, there is a coexistence of an elec- 
trophoretic and a dielectrophoretic force in the presence of 
a nonuniform AC electric field. The DEP force always points 
toward the region of high field gradient. It does not oscil- 
late with the change of direction of the field. In contrast, the 
electrophoretic force points along the direction of field, and 
hence is oscillatory under the same conditions. How to sepa- 
rate the DEP force from the electrophoretic force is a question 
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FIG. 8: Dispersion strength versus the characteristic frequency for 
different medium conductivities. Parameters: £ = 0.12V, r\ = 1.0 x 
10~ 3 Kg/(ms), R/D = 1.03, ei = 2.25e , u = 0.033 C/m 2 , E = 
0.014 Sm 2 /mol. The lines are drawn as a guide to the eye. 

of interest in many experimental setups [^| ^]. In different 
frequency ranges, either the electrophoretic force or the DEP 
force dominates, and the transition from one to the other oc- 



curs at a frequency f tr , which has been approximately deter- 
mined Here, we have chosen a frequency region where 
electrophoretic effects are negligible and the DEP force dom- 
inates. In addition, although we are at finite temperature, 
Brownian motion is not included in our analysis. In experi- 
ments Brownian motion is always present and has posed diffi- 
culties in dielectrophoresis of submicrometer particles. How- 
ever, with current techniques it is possible to access also this 
range [gg. 

One of the interesting questions is what happens, when the 
volume fraction of the suspension becomes large. It turns out 
that it is possible to extend our approach by taking into ac- 
count local field effects which may modify the DEP crossover 
spectrum. Work is in progress to address these questions. In 
addition to dielectrophoresis, the extension of the present ap- 
proach is also of interest from the point of view of electroro- 
tation. 

To summarize, using the multiple image method, we have 
been able to capture mutual polarization effects of two ap- 
proaching particles in an electrolyte. Using spectral represen- 
tation theory, we derived an analytic expression for the DEP 
force, and using that the crossover frequency was determined. 
From the theoretical analysis, we find that the mutual polariza- 
tion effects can change the crossover frequency substantially. 
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